#include "estimador.h"

int estimador_varianza( const float *m , const SIZE_T n, const SIZE_T col , float *result ) {
	const float *upper_limit;				// t0
	upper_limit = m + n*n;				// load

	float total = 0;						// f20
	float total_de_cuadrados = 0;

	m += col;	
	
	for (;  m < upper_limit; m+=n) {
		
		float f4 = *m;
		
		total += f4;					// load
		total_de_cuadrados += f4 * f4;
	};
	
	*result = total_de_cuadrados / n - total*total / n / n ;

	return 0;
}




